Technical  Report  HL-92-7 
September  1995 


WwW 


Min 


US  Army  Corps 
of  Engineers 

Waterways  Experiment 
Station 


Houston-Galveston  Navigation  Channels, 
Texas  Project 


Report  4 

Three-Dimensional  Numerical  Modeling 
of  Hydrodynamics  and  Salinity 


by  R.  C.  Berger,  Robert  T.  McAdory,  Joseph  H.  Set 
William  D.  Martin,  WES 
Larry  H.  Hauck,  Jones  and  Neuse  Engineering 


k  DtG  0  8 


Approved  For  Public  Release;  Distribution  Is  Unlimited 


lOTC  Qtf AJuTH  INSPBCrEl^  3 


Prepared  for  U.S.  Army  Engineer  District,  Galveston 


19951204  101 


Technical  Report  HL-92-7 
September  1995 


Houston-Galveston  Navigation  Channeis, 
Texas  Project 

Report  4 

Three-Dimensional  Numerical  Modeling 
of  Hydrodynamics  and  Salinity 

by  R.  C.  Berger,  Robert  T.  McAdory,  Joseph  H.  Schmidt,  William  D.  Martin 

U.S.  Army  Corps  of  Engineers 
Waten/vays  Experiment  Station 
3909  Halls  Ferry  Road 
Vicksburg,  MS  39180-6199 

Larry  H.  Hauck 

Jones  and  Neuse  Engineering 
Austin,  TX  77553 


Accesion  For 

NTIS  CRA&I  to 

OTIC  TAB  □ 

Unannounced  □ 

Justification 

By 

Distfib 

jtion  / 

Availability  Codes 

Dist 

Avail  and/or 

Special 

Report  4  of  a  series 

Approved  for  public  release;  distribution  is  unlimited 


Prepared  for  U.S.  Army  Engineer  District,  Galveston 
Galveston,  TX  77550 


N 


US  Army  Corps 
of  Engineers 

Waterways  Experiment 
Station 


HEADQUARTERS 


PDR  WFOrtJATlON  OCJWTACT  : 

PUBLIC  AF=FA1RS  OFFICE 
U.  S.  ARMY  aJGINEER 
WATERWAYS  EXPERJMEKT  STATION 
3909  HALLS  FERRY  ROAD 
VICKSBURG,  Mississippi  39160-6199 
PHONE :  {601)63A-2502 


AREA  OF  R^RVATTCN  .  2.7  *q  Vn 


Waterways  Experiment  Station  Cata!oging-in-Pub!ication  Data 

Houston-Galveston  Navigation  Channels,  Texas  Project.  Report  4, 
Three-dimensional  numerical  modeling  of  hydrodynamics  and  salinity/ 
by  R.C.  Berger ...  [et  al.] ;  prepared  for  U.S.  Army  Engineer  District, 
Galveston. 

482  p. :  ill. :  28  cm.  —  (Technical  report ;  HL-92-7  rept.4) 

Includes  bibliographical  references. 

Report  4  of  a  series. 

1 .  Salinity  —  Texas  —  Galveston  Bay.  2.  Galveston  Bay  (Tex.) 

3.  Houston  Ship  Channel  (Tex.)  I.  Berger,  Rutherford  C.  II.  United 
States.  Army.  Corps  of  Engineers.  Galveston  District.  III.  U.S.  Army 
Engineer  Waterways  Experiment  Station.  IV.  Hydraulics  Laboratory 
(U.S.  Army  Engineer  Waterways  Experiment  Station)  V.  Title:  Three- 
dimensional  numerical  modeling  of  hydrodynamics  and  salinity. 

VI.  Series:  Technical  report  (U.S.  Army  Engineer  Waterways  Experiment 
Station) ;  HL-92-7  rept.4. 

TA7  W34  no.HL-92-7  rept.4 


Contents 


Preface  .  v 

Conversion  Factors,  Non-SI  to  SI  Units  of  Measurement . vi 

1 —  ^Introduction .  1 

Background  .  1 

Objectives  .  2 

2 —  ^Development  of  the  Model  .  4 

Development  of  the  Scope  of  Work .  4 

The  Model .  5 

Data  Requirements  .  8 

Verification  . 11 

3 —  ^Testing  and  Results . 15 

Introduction  . 15 

Boundary  Conditions . 16 

Test  Procedures . 19 

Results . 19 

Time-series  . 19 

Isohaline  plots  . 24 

Mitigation  plan . 25 

4 —  Summary  and  Conclusions . 26 

References  . 28 

Plates  1-425 

Appendix  A:  The  Hydrodynamic  Code . A1 

Appendix  B:  Tabulation  of  Freshwater  Inflows  . B1 

SF  298 

DXIC  QUaXaII*^  3 

List  of  Figures 


Figure  1.  Location  map  of  Houston-Galveston  Navigation 
Channel  Project . 


3 


Figure  2.  Existing  conditions . 7 

Figure  3.  Phase  I  conditions . 8 

Figure  4.  Phase  II  conditions . 9 

Figure  5.  NED  conditions  .  10 

Figure  6.  Mitigation  conditions  .  11 

Figure  7.  2>D  and  3-D  delineation  .  12 

Figure  8.  Basic  long-term  survey  data  equipment  location .  13 

Figure  9.  Supplementary  long-term  collection  equipment 

location  map  .  14 

Figure  10.  Total  freshwater  inflows .  18 

Figure  11.  Time-series  points  . 20 


Preface 


The  Ihree-dimensional  numerical  modeling  testing  of  hydrodynamic  and 
salinity  conditions  for  the  Houston-Galveston  Navigation  Channels,  Texas 
Project,  as  documented  in  this  report,  was  performed  for  the  U.S.  Army  Engi¬ 
neer  District,  Galveston  (SWG),  by  personnel  of  the  Hydraulics  Laboratory 
(HL)  of  the  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES),  The 
testing  program  reported  herein  occurred  from  February  1992  through  January 
1995.  The  study  began  in  October  1990.  The  Houston-Galveston  Navigation 
Channels,  Texas  Project,  study  was  authorized  by  House  Document  350, 

85th  Congress,  2nd  Session,  and  by  House  Document  427,  86th  Congress, 

2nd  Session.  The  subject  study  was  authorized  by  SWG  in  conjunction  with 
the  study  cost-sharing  partner,  the  Ports  of  Houston  and  Galveston,  Texas,  on 
3  April  1990. 

This  is  Report  4  of  a  series.  Report  1  describes  the  data  collection. 

Report  2  is  the  two-dimensional  modeling  of  hydrodynamics  for  the  navigation 
study,  and  Report  3  describes  the  verification  of  the  three-dimensional  model 
for  hydrodynamics  and  salinity. 

The  study  was  conducted  under  the  general  supervision  of  Messrs.  F.  A. 
Herrmann,  Jr.,  Director,  HL;  R.  A.  Sager,  Assistant  Director,  HL;  W.  H. 
McAnally,  Jr.,  Chief,  Waterways  and  Estuaries  Division  (WED),  HL;  and 
W.  D.  Martin,  Acting  Chief,  Hydro-Science  Division  (HSD),  HL. 

This  work  was  performed  and  the  report  prepared  by  Dr.  R.  C.  Berger  and 
Dr.  R.  T.  McAdory,  WED;  Messrs.  J.  H.  Schmidt  and  Martin,  HSD;  and 
Mr.  L.  H.  Hauck,  previously  associated  with  WED  and  Jones  and  Neuse  Engi¬ 
neering.  Much  of  the  daily  activities  such  as  checking  run  status,  plotting,  and 
data  storage  was  performed  by  Messrs.  Jay  Hardy  and  Lenwaski  Campbell, 
WED,  and  Ms.  Cassandra  Gaines,  HSD. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was 
Dr.  Robert  W.  Whalin.  Commander  was  COL  Bruce  K.  Howard,  EN. 


The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication, 
or  promotional  purposes.  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 


Conversion  Factors, 
Non-SI  to  SI  Units  of 
Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units 
as  follows: 


Multiply 

By 

To  Obtain 

acre -feet 

1,233.489 

cubic  meters 

acres 

4,046.873 

square  meters 

cubic  feet 

0,02831685 

cubic  meters 

feet 

0.3048 

meters 

miles  (U.S,  nautical) 

1.852 

kilometers 

miles  (U.S.  statute) 

1.609347 

kilometers 

1  Introduction 


Background 

In  the  19lh  century  Galveston  Bay  existed  as  an  immense  shallow  estuary 
accessible  to  light-draft  vessels.  Conditions  then,  as  today,  created  a  fertile 
source  for  shell  and  fin  fish.  Areas  of  the  bay,  such  as  Red  Fish  Reef,  were, 
according  to  verbal  records,  so  shallow  that  during  low  tides  cattle  could  be 
driven  across  the  bay.  Subsequently,  as  the  Texas  Gulf  Coast  developed, 
channels  were  constructed  to  provide  access  to  Galveston,  Texas  City,  and 
Houston,  TX,  for  increasingly  deeper  draft  vessels.  A  brief  chronology  of  the 
deepening  is  shown  in  Table  1. 


Table  1 

Chronology  of  Houston  Channel  Deepening,  Galveston  Bay  (from 
Alperin  1977) 

Date 

Channel  Minimum  Depth,  ft 

Channel  Width,  ft 

1851 

4 

- 

1870 

4 

70 

1874 

9 

120 

1889 

12 

100 

1893 

14 

100 

1903 

18.5 

150 

1914 

25 

150 

1926 

30 

250 

1932 

32 

250 

1936 

34 

400 

1948 

36 

400 

i 

1964 

40 

400 

The  currently  proposed  Houston-Galveston  Navigation  Project  is  the  latest 
attempt  to  maintain  the  competitiveness  of  this  area  in  the  world  trade  market. 
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Houston  is  currently  the  number  two  port  in  the  United  States  based  on  ton¬ 
nage  handled,  rated  behind  only  New  Orleans  (’The  U.S.  Waterway  System 
Fact  Card"  1991). 

The  project  was  originally  proposed  for  construction  in  two  construction 
phases.  In  Phase  I,  the  Houston  Ship  Channel  (HSC)  will  be  enlarged  to  a 
depth  of  45  ft^  and  widened  to  a  width  of  530  ft,  and  the  Galveston  Channel 
to  a  depth  of  45  ft  and  a  width  of  450  ft.  In  Phase  II,  both  channels  were  to 
be  deepened  to  50  ft  and  HSC  widened  to  600  ft.  Current  plans  do  not 
include  the  Phase  II  enlargement.  Both  channel  depth  conditions  were 
analyzed  before  this  decision  was  finalized  and  are  thus  included  in  this  report. 
In  the  remainder  of  the  report,  the  45-ft-deep  channel  is  designated  as  Phase  I 
and  the  50-ft-deep  channel  as  Phase  II.  The  study  area  is  shown  in  Figure  1. 


Objectives 

There  are  several  related  objectives  of  this  study; 

a.  Develop  an  accurate,  verified  three-dimensional  hydrodynamic  and 
salinity  model  of  Galveston  Bay  for  the  conditions  existing  in  late  1990. 

b.  Using  this  model,  evaluate  the  effects  of  changes  due  to  proposed 
channel  deepening  on  circulation  and  salinity  within  the  bay. 

c.  Evaluate  the  effects  of  the  channel  deepening  under  several  predicted 
future  changes  in  freshwater  inflow. 

d.  Provide  the  output  from  these  model  runs  as  input  to  a  comprehensive 
model  for  evaluating  oyster  population  response  to  the  changed 
conditions  modeled. 


^  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  units  is  found  on  page  vi. 


2 


Chapter  1 


Introduction 


Chapter  1  Introduction 


3 


Figure  1.  Location  map  of  Houston-Galveston  Navigation  Channel  Project 


2  Development  of  the  Model 


Development  of  the  Scope  of  Work 

From  the  inception  of  the  studies  at  the  U.S.  Army  Engineer  Waterways 
Experiment  Station  (WES),  interagency  coordination  was  an  integral  part  in 
developing  the  scope  of  work  and  reporting  progress.  The  Interagency 
Coordination  Team  (ICT)  provided  the  forum  for  discussions  of  the  scope  of 
work  for  the  WES  studies  in  addition  to  the  many  other  topics  and  studies  that 
are  part  of  the  Houston-Galveston  Navigation  Project.  Member  organizations 
of  the  ICT  include  U.S.  Army  Engineer  District,  Galveston;  Houston  Port 
Authority;  U.S.  Fish  and  Wildlife  Service;  U.S.  Environmental  Protection 
Agency;  National  Marine  Fisheries  Service;  Texas  Governor’s  Office;  Texas 
Parks  and  Wildlife  Department;  Galveston  Bay  National  Estuary  Program;  and 
Galveston  Wharves.  Over  the  period  January  1990  through  September  1990, 
WES  and  ICT  went  through  an  interactive  process,  which  culminated  in  a  final 
scope  of  work  from  WES  dated  September  27,  1990.  In  the  process  of 
developing  the  scope  of  work,  WES  personnel  attended  the  ICT  meetings  on 
approximately  a  monthly  basis  for  presentation  of  the  evolving  scope,  and  also 
to  present  progress  on  studies  that  had  prior  endorsement  from  the  ICT,  for 
example,  the  velocity  development  for  the  ship  simulator  studies  and  initial 
phases  of  the  field  investigations. 

As  part  of  the  interagency  coordination,  the  WES  staff  provided  a  presen- 
Lation  on  April  12,  1990,  to  the  Ocean  and  Estuarine  Physics  Branch,  National 
Ocean  Service  (NOS),  a  sister  organization  to  the  National  Marine  Fisheries 
Service,  both  in  the  National  Oceanic  and  Atmospheric  Administration 
(NOAA).  In  this  presentation,  a  comprehensive  description  of  the  state-of-the- 
art  three-dimensional  (3-D)  model,  RMAIO-WES,  was  provided  to  the  NOS 
staff.  A  series  of  phone  conversations  and  letters  ensued  with  NOS  agreeing 
that  the  RMAIO-WES  model  should  be  capable  of  addressing  the  objectives  of 
the  Galveston  Bay  Study. 

Final  agreement  to  the  field  investigation  plan  was  reached  during  a  meet¬ 
ing  on  September  13,  1990,  of  U.S.  Army  Corps  of  Engineers,  NOAA  (NOS 
and  National  Marine  Fisheries),  and  Council  on  Environmental  Quality,  which 
resulted  in  a  statement  of  agreement.  The  agreement  to  the  field  investigation 
involved  the  following  two  major  items:  (a)  extension  of  the  long-term 
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instrument  deployment  from  120  days  to  180  days  (July  19,  1990,  to  mid- 
January,  1991),  and  (b)  redeployment  of  several  of  the  moored  current  meters 
and  salinometers  to  new  locations  in  October  or  November.  These  field 
investigation  modifications  allowed  extension  of  instrument  deployment  into 
the  winter  to  record  frontal  passages,  which  are  important  meteorological 
events  influencing  bay  circulation.  Also,  the  instrument  redeployment  allowed 
dual  emphasis  of  data  gathering — the  HSC  for  the  first  half  of  the  study  and 
the  shallows,  especially  Trinity  Bay,  for  the  last  half.  The  results  of  the  field 
data  collection  effort  are  presented  in  Fagerburg  et  al.  (1994). 

Concurrent  with  the  development  of  the  field  investigation  scope  of  work 
was  the  development  of  the  3-D  (RMAIO-WES)  modeling  scope  of  work.  The 
number  and  length  of  the  RMAIO-WES  applications  were  largely  determined 
by  the  need  to  simulate  the  same  hydrologic  conditions  as  in  the 
Environmental  Impact  Statement  (EIS)  and  to  provide  sufficient  hydraulic  data 
to  assess  project-induced  impacts  on  oysters. 

By  the  summer  of  1990,  ICT  had  agreed  to  the  scope  of  work  for  a  numer¬ 
ical  oyster  model  to  evaluate  project-induced  impacts  to  oysters.  Oysters  were 
selected  as  the  single  species  whose  fate  is  most  indicative  of  the  general  bay 
condition  and  its  impact  on  other  important  species.  The  oyster  model  studies 
were  directed  by  Dr.  Eric  Powell  of  Texas  A&M  University.  Because  the 
oyster  model  used  currents  and  salinities  from  RMAIO-WES,  the  WES  3-D 
model  studies  were  strongly  linked  to  the  oyster  model  studies.  The  duration 
of  several  key  RMAIO-WES  simulations  was  dictated  largely  by  the  oyster 
model  needs  for  complete  information  during  the  oyster’s  life  cycle.  Conse¬ 
quently,  for  these  key  hydrologic  sequences,  simulations  of  12-month  duration 
were  conducted.  Other  simulations  were  run  for  the  period  January-September, 
which  encompasses  the  oyster’s  annual  active  stage. 


The  Model 


The  task  of  modeling  Galveston  Bay  in  three  dimensions  while  varying  the 
boundary  conditions  and  geometries  required  use  of  a  sophisticated  numerical 
model.  Several  models  were  considered  before  the  RMAIO-WES  code  was 
chosen.  This  finite  element  code  allowed  easy  modification  of  the  geometry 
and  bathymetry.  The  challenge  lay  in  the  length  of  the  simulations.  The 
model  had  previously  been  used  to  model  real-time  durations  of  several  days; 
yet  the  3-D  study  required  several  runs  of  1  year  and  others  of  9  months. 

The  unstructured  mesh  capability  of  RMAIO-WES  allowed  flexible  resolu¬ 
tion  that  matched  small-  and  large-scale  features.  Several  disposal  sites  for 
dredged  material  are  included  for  both  phases.  These  sites  vary  from  as  little 
as  100  acres  to  as  much  as  1,500  acres.  They  are  shown  on  the  computational 
grids  for  Phases  I  and  II,  to  be  discussed  later  in  this  section.  The  National 
Economic  Development  (NED)  plan  has  the  same  channel  configuration  as 
Phase  I  but  with  submerged  disposal  material  placement  instead  of  the  islands 
in  Phase  I. 
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The  RMAIO-WES  code  is  a  Galerkin-bascd  finite  element  solution  to 
simulate  3-D  unsteady  opcn-channcl  How.  The  model  was  originally 
developed  by  Dr.  Ian  King  of  Resource  Management  Associates  (King  1988) 
and  extensively  modified  by  the  WES  staff.  The  code  represents  3-D  hydro¬ 
dynamics  using  conservation  of  fluid  mass,  horizontal  momentum,  and 
salinity/temperature  transport  equations.  As  is  typical  of  shallow-water 
models,  the  vertical  accelerations  are  assumed  to  be  negligible  (the  hydrostatic 
assumption)  and  so  the  vertical  velocity  is  calculated  through  mass  conser¬ 
vation.  In  the  interest  of  computational  efficiency  the  code  also  simulates  one¬ 
dimensional  (1-D)  and  two-dimensional  (2-D)  flow  as  well  as  transitions 
between  1-,  2-  and  3-D.  The  code  is  implicit  in  time  and  resolves  the  non- 
linearities  via  Newton-Raph.son  iteration. 

In  a  3-D  model,  or  a  laterally  integrated  2-D  model,  the  vertical  extent  of 
the  domain  is  not  known  until  the  depth  is  calculated,  and,  of  course,  one 
cannot  make  the  ctilculations  until  the  computational  grid  is  developed.  This 
apparent  impasse  is  avoided  by  transforming  the  domain  at  each  time-step  to  a 
static  grid.  The  particular  transformation  used  maps  the  water  surface  to  a 
constant  elevation,  the  bed  is  unchanged,  and  all  elevations  in  between  are 
stretched  proportionally.  Five  such  grids  were  used  for  this  study.  These 
represented  the  1990,  or  existing,  conditions  (Figure  2);  the  45-ft-dcep  channel 
(Phase  I)  with  17  intermediate  mitigation  islands  added  (Figure  3);  the  50-ft- 
deep  channel  (Phase  II)  with  intermediate  mitigation  islands  added  (Figure  4); 
the  NED  plan  (Figure  5);  and  the  final  mitigation  plan  with  45-ft-dccp 
channels  (Figure  6).  The  characteristics  of  the  grids  are  summarized  in  the 
following  tabulation.  The  mitigation  plan  run  was  made  lo  supply  the  oyster 
model  with  data  representative  ot  the  selected  plan  of  improvement.  The 
oyster  model  would  use  these  data  to  determine  the  size  and  location  of  oyster 
reefs  needed  to  mitigate  the  adverse  impacts  of  the  selected  plan. 


Condition  Modeled 

Surface 

Nodes 

Surface 

Elements 

Total  Nodes 

Total  Elements 

Existing 

6.192 

1,996 

12,270 

5,112 

45-ft  channel 

6,902 

2,278 

14,477 

6,250 

50-ft  channel 

6,902 

2,278 

14,477 

6,250 

NED 

7,027 

2,333 

14,990 

1 

'6,421 

Mitigation 

7,414 

2,451 

16,470 

7,045 

The  model  is  configured  as  largely  3-D  with  the  exception  of  the  West  Bay 
west  of  the  causeway.  The  3-D  portion  consists  of  one  element  (three  nodes) 
in  the  shallow  areas  of  the  bay  and  up  to  four  elements  (nine  nodes)  in  the 
channel.  Figure  7  depicts  the  2-D  and  3-D  portions  of  the  grid.  These  por¬ 
tions  are  the  same  for  each  case  except  where  mitigation  islands  are  present  as 
solid  land  mass. 

The  salinity/density  relationship  in  the  model  is  based  on  Pritchard  (1982). 
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Figure  2.  Existing  conditions 


The  vertical  turbulence  is  a  combination  of  the  Mellor-Yamada  Level  II 
(Mellon  and  Yamada  1982;  Adams  and  Weatherly  1981)  and  Hendcrson-Sellers 
(1984).  The  wind  stress  is  a  quadratic  relationship  using  the  specific  coeffi¬ 
cients  developed  by  Wu  (1980).  The  effect  of  wind  waves  on  vertical  mixing 
is  simulated  in  the  manner  of  the  numerical  model  CE-QUAL-W2  (U.S.  Army 
Engineer  Waterways  Experiment  Station,  1986).  Details  of  the  model  are 
presented  in  Appendix  A. 
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Figure  3.  Phase  I  conditions 


Data  Requirements 

Any  model  requires  extensive  field  data  to  ensure  that  the  model  is 
accurately  reproducing  conditions  observed  in  the  field.  To  this  end,  tides, 
salinities,  and  current  magnitudes  and  directions  at  several  depths  were  col¬ 
lected  at  28  different  locations  during  the  period  July  1990  through  January 
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Figure  4.  Phase  II  conditions 


1991  (Figure  8).  These  data  represent  3,105  days  of  recorded  data  during  this 
time  frame.  Additionally,  wind  speed  and  direction  were  recorded  for  this 
period  at  sta  SlO.l  (Figure  9)  and  at  Houston  Intercontinental  Airport.  A  full 
description  of  these  data  is  presented  in  Fagerburg  et  al.  (1994).  Use  of  these 
data  is  discussed  in  Berger  et  al.  (in  preparation). 

Geometric  and  bathymetric  data  for  the  bay  were  obtained  from  NOAA 
Charts  numbered  11,323,  11,326,  and  11,332.  Additionally,  existing  channel 
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Figure  5.  NED  conditions 


conditions  were  provided  by  the  U.S.  Army  Engineer  District,  Galveston.  The 
District  also  provided  channel  alignments,  depths,  and  side  slopes  for  the  two 
deepening  phases,  some  of  which  had  been  determined  from  the  ship  simulator 
covered  in  Hewlett  (^1994).  Freshw'ater  inflows  lor  the  various  scenarios 
were  developed  jointly  by  the  Galveston  District  and  the  Texas  Water 
Development  Board.  These  flows  considered  future  water  use,  source  of  this 
water,  known  projects  such  as  the  Wallisvillc  Dam  on  the  Trinity  River,  and 
projected  future  interbasin  water  diversions. 
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Figure  6.  Mitigation  conditions 


Verification 


The  period  of  data  collection  was  from  19  July  1990  to  15  January  1991. 
Data  were  collected  by  moored  salinity  and  velocity  meters,  as  well  as  tide 
gauges  mounted  throughout  the  system,  during  the  bulk  of  this  lime.  Within 
this  6-month  period  a  short-duration  intense  field  collection  was  conducted  on 
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Figure  7.  2-D  and  3-D  delineation 


19-20  July  1990.  These  data  included  velocity  and  salinity  measurements  over 
depth  at  five  ranges  along  the  navigation  channel. 

The  basic  philosophy  adopted  in  the  verification  process  was  to  make  any 
appropriate  adjustments  to  the  model  to  reasonably  match  the  short  intense 
survey  period,  then  make  the  long-term  simulations  unaltered  and  compare 
them  to  measured  results  as  a  check  on  the  model's  performance.  This  effort 
is  fully  reported  in  Berger  et  al.  (in  preparation). 
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Figure  8.  Basic  long-term  survey  data  equipment  locations  (stations  S1 .0,  S2.0,  S5.5,  S6.0, 
S6.1,  S8.1,  S9.0,  S11.1,  SI  2.0,  and  S15.1  record  salinity  as  well) 
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3  Testing  and  Results 


Introduction 

This  chapter  describes  the  basic  test  conditions,  the  manner  in  which  the 
model  was  run,  results,  and  discussion  of  these  results.  The  basic  test  condi¬ 
tions  consist  of  duration  of  simulation,  the  bathymetry/geometry,  hydrology, 
and  the  tidal  and  salinity  boundary. 

A  total  of  60  possible  scenarios  was  considered  tor  RMAIO-WES  simula¬ 
tion.  Five  geometric  conditions  were  evaluated,  and  each  could  be  evaluated 
for  four  hydrologic  conditions.  These  hydrologic  conditions  are  each  further 
subdivided  into  three  representative  flow  levels.  The  five  geometric  conditions 
are  as  follows: 

a.  Existing  channel. 

b.  National  Economic  Development  Plan  (NED)  channel  and  dredged 
material  disposal  plan  (the  channel  condition  in  the  EIS). 

c.  Phase  I  channel  improvements  (45-ft  depth)  with  intermediate  dredged 
material  beneficial  use  plan. 

d.  Phase  II  channel  improvements  (50-ft  depth)  with  intermediate  dredged 
material  beneficial  use  plan. 

e.  Final  mitigation  (45-ft-depth)  channel  with  final  dredged  material 
beneficial  use  plan. 

The  following  four  hydrologic  conditions  reflect  projected  future  hydrologic 
and  water  use  patterns  considering  planned  reservoirs,  water  importation, 
groundwater  usage  changes,  and  population  growth: 

a.  Present  (year  1990). 

b.  Year  1999  (with  proposed  Wallisville  saltwater  barrier). 

c.  Year  2024. 
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d.  Year  2049. 


For  all  four  hydrologic  conditions,  three  time-varying  freshwater  inflow 
levels  representing  an  annual  sequence  of  high,  median,  or  low  inflow  were 
considered.  These  flows  are  tabulated  by  location  and  month  as  acre-feet  in 
Appendix  B.  A  subset  of  all  these  possible  scenarios  was  selected  to  produce 
a  manageable  analysis  that  would  meet  the  objectives  of  this  study. 

The  Gulf  of  Mexico  salinity  varies  by  season.  Mississippi  River  outflow 
hugs  the  Texas  coast  and  lowers  salinities  during  certain  seasons.  The  work  of 
Cochrane  and  Kelly  (1986)  was  consulted  and  15-year  monthly  averages  for 
Gulf  salinity  obtained.  These  were  used  for  the  Gulf  boundary  salinities  for 
each  month  run  in  the  model.  These  are  shown  in  the  following  tabulation: 


Month 

Average  Monthly 

1957-1971 

Salinity,  ppt 

January 

29.1 

February 

29.1 

March 

28.4 

April 

26.2 

May 

24.9 

June 

28.1 

July 

31.6 

August 

33,3 

September 

28.7 

October 

28.5 

November 

29.5 

December 

29.4 

Boundary  Conditions 

The  series  ot  tests  simulated  in  the  model  arc  shown  in  the  following 
tabulation. 

These  runs  produced  salinity  values  over  the  9-month  period  for  a  select 
number  of  locations,  based  upon  the  hydrodynamic  model  runs  at  other  flow 
hydrologies.  This  approach  is  covered  in  Hsich  (in  preparation),  and  the 
results  of  those  runs  presented  in  that  report. 

The  period  of  primary  interest  for  oysters  and  other  species  is  the  summer 
months.  Thus,  the  Y-runs  were  considered  to  be  long  enough  to  provide  data 
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Geometry/Bathymetry 

Hydrology 

Present 

1999 

2024 

2049 

L 

M 

H 

L 

M 

H 

L 

M 

H 

L 

M 

H 

Existing 

X 

X 

X 

Y 

S 

- 

Y 

Y 

Y 

Y 

Y 

Y 

Phase  1 

X 

X 

X 

Y 

S 

- 

Y 

Y 

Y 

Y 

Y 

Y 

Phase  li 

- 

- 

- 

- 

- 

- 

Y 

Y 

Y 

Y 

Y 

Y 

NED 

- 

- 

- 

- 

- 

- 

Y 

- 

- 

- 

- 

- 

Mitigation 

- 

- 

- 

- 

- 

- 

- 

X 

- 

- 

- 

- 

Note:  X  =  runs  1  year  long,  January-December. 

Y  =  runs  9  months  long,  January-September. 

S  =  runs  generated  using  systems  model  approach. 

for  oyster  model  results  tor  9-monlh  runs  and  12  months  of  data  from  the 
X-runs. 


The  various  hydrologic  eras  represent  predicted  water  use  and  hydrologic 
patterns  expected  over  the  life  of  the  project.  These  are  described  in 
U.S.  Army  Engineer  District,  Galveston  (1987),  Texas  Water  Development 
Board  (1990),  and  Metcalf  and  Eddy  with  Ekistics  Corporation  (1989).  The 
1999  era  includes  the  Wallisville  saltwater  barrier  project  on  the  Trinity  River 
and  the  redistribution  of  Hows  accompanying  the  release  operations.  The 
2024  era  represents  a  time  in  which  the  projected  growth  of  the  Houston 
metropolitan  area  forces  a  diversion  of  water  from  the  Trinity  River  to  the  city 
and  ultimately  into  Galveston  Bay.  This  is  generally  regarded  as  the  era  that 
represents  a  worst  case  for  the  bay  system.  The  2049  era  is  at  the  end  of 
estimated  project  life.  With  the  continued  projected  growth  of  the  Houston 
area,  additional  flow  is  diverted  into  the  Galveston  Bay  system  trom  rivers  or 
water  sources  outside  the  Galveston  Bay  watershed;  i.e.,  it  is  projected  that 
there  will  be  more  fresh  water  flowing  into  the  Galveston  system  than  in  2024. 
The  low,  medium  and  high  are  typical  flow  years  of  that  nature  in  the  pro¬ 
jected  year  2049.  These  inflows  were  synthesized  by  the  Texas  Water  Devel¬ 
opment  Board  and  Galveston  District  to  reflect  projected  typical  low-, 
medium-,  and  high-flow  years  in  the  various  eras.  The  data  are  given  as 
monthly  averages  and  depicted  in  Figure  10. 

Two  Houston  Power  and  Light  Company  electric  power  generation  plants 
were  included  in  the  model:  the  Cedar  Bayou  and  Robinson  plants.  Monthly 
power  plant  intake  rates  for  the  year  1990  were  used  for  all  testing.  The 
model  used  the  salinity  values  obtained  at  the  intake  as  the  specified  boundary 
condition  at  the  discharge. 

Flows  from  Buffalo  Bayou,  San  Jacinto  River,  and  Trinity  River  were 
enforced  as  specified  discharge  boundary  conditions  across  particular  boundary 
nodes  with  a  salinity  value  of  zero  concentration.  On  the  small  streams. 
Chocolate,  Highland,  Oyster,  Robinson,  Cedar,  Dickinson,  and  Live  Oak 
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Figure  10.  Total  freshwater  inflov/s 

Bayous  and  Clear  Creek,  the  flow  rates  arc  so  small,  on  the  order  of  200  cl's, 
that  if  one  implemented  a  zero  salinity  boundary  condition  specified  at  these 
boundary  nodes,  the  fliiid  nearby  would  not  be  significantly  diluted.  An 
extraordinary  amount  of  resolution  would  then  be  required  to  resolve  these 
gradients.  Instead  the  How  was  added  to  appropriate  elements  in  the  same 
manner  as  one  would  input  rain,  for  example.  This  allowed  the  area  to  be 
diluted  without  requiring  excessive  resolution. 

The  freshwater  inllows,  power  plant  intake/discharge,  and  salinity  boundary 
data  were  updated  every  day  in  the  simulation.  To  prevent  extremely  rapid 
variation  in  these  values,  a  9-day  moving  average  wais  used  to  slightly  smooth 
the  data.  This  was  particularly  important  in  the  freshwater  inllows  since  they 
were  given  only  as  monthly  values  and  changes  between  months  could  have 
been  large. 

The  tidal  boundary  data  were  from  NOS.'NOAA  station  877-1510  at  Plea¬ 
sure  Pier  for  calendar  year  1984.  This  year  was  found  to  be  a  typical 
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temperature  year,  which  was  felt  to  be  appropriate  for  the  oyster  modeling 
work.  The  model  boundary,  which  is  offshore  about  26  miles,  used  these  data 
shifted  forward  in  time  by  1.31  hr.  This  matched  the  phase  at  Pleasure  Pier. 
Periods  of  less  than  3  hr  were  removed  by  filtering,  preventing  noise  in  the 
results. 

The  1984  wind  data  were  obtained  from  the  National  Weather  Service  for 
Houston  Intercontinental  Airport.  From  a  correlation  with  data  collected  at 
sta  SlO.l  in  Galveston  Bay  during  the  verification  period,  these  wind  data 
were  transposed  to  the  bay  system. 


Test  Procedures 


The  original  boundary  condition  information  was  produced  on  a  daily  basis 
for  all  freshwater  inflows  and  power  plants.  The  wind  magnitude  and  direc¬ 
tion,  tide  elevation,  and  salinity  were  hourly.  The  largest  time-step  taken  was 
1  hr.  If  the  convergence  suggested  that  smaller  time-steps  should  be  run,  the 
boundary  data  were  interpolated  to  either  half-  or  quarter-hours,  as  necessary, 
and  half-  or  quarter-hour  time-steps  were  run.  The  problem  solved  was  non¬ 
linear  in  that  several  of  the  unknown  variables  were  written  as  products.  Thus 
a  system  to  linearize  the  problem  and  iterate  each  time-step  to  final  solution 
was  required.  The  method  chosen  was  the  Newton-Raphson  method.  The 
number  of  iterations  in  the  Newton-Raphson  scheme  was  either  3  or  4,  as 
needed.  Most  of  these  procedures  were  automatically  implemented  by  the 
model  or  the  scripts  (computer-run  control  programs)  that  ran  the  model. 

Hot-start  files  were  generated  regularly,  which  consisted  of  all  values  at  the 
computational  nodes  at  a  point  in  time.  This  file  could  then  be  run  to  restart 
the  model  “hot,”  which  requires  no  spin-up  time  and  basically  picks  up  where 
the  previous  run  ended.  A  file  of  select  nodal  data  at  all  depths  was  saved  and 
plotted  for  quality  control,  and  all  hydrodynamic  and  salinity  data  were  saved 
on  file  for  eventual  export  for  use  in  the  oyster  model. 


Results 


Time-series 

The  time-series  results  for  Phase  I  versus  existing  bathymetry  (present 
hydrology,  1990)  flow  are  shown  in  Plates  1-48.  The  bottom  salinity  results 
for  low  flow  are  in  Plates  1-16,  medium  flow  in  Plates  17-32,  and  high  flow  in 
Plates  33-48.  The  overall  trends  are  as  noted  in  the  section,  "Isohaline  plots," 
with  the  most  significant  salinity  changes  occurring  in  the  upper  bay  west  of 
the  channel.  The  time-series  results  yield  information  about  the  time  distri¬ 
bution  and  duration  of  the  salinities  at  particular  locations,  shown  in  Figure  11. 
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TRINITY  RIVER 


Figure  11.  Time-series  points 


Table  2  summarizes  all  of  the  time-history  data  grouped  into  2-month  seg¬ 
ments  (April-May,  June-July,  August-September,  October-November).  The 
values  given  are  in  terms  of  the  average  salinity  increase:  a  value  of  -1,  for 
example,  indicates  the  plan  showed  an  average  decrease  in  salinity  between  0 
and  1  ppt.  A  value  of  0  means  increases  between  0  and  1  ppt,  a  value  of 
1  indicates  the  increase  is  between  1  and  2  ppt,  etc.  Essentially,  this  amounts 
to  rounding  down  to  the  nearest  integer  value  of  salinity  change.  The  periods 
of  less  than  14  days  (350  hr)  have  been  removed  by  filtering.  As  a  result  of 
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Table  2 

Summary  of  Time-History  of  Salinity  Change  Results,  2-Month  Averages, 
Bottom  Elevation 


Salinity,  ppt 


Trinity 

Hydrology  Months  Point  2 


1990  Low 


Channel 

East  Bay 

South 

West 

West  Bay 

Point  3 

Point  7 

Point  8 

Point  12 

Phase  I 


Chapter  3  Testing  and  Results 


21 


Table  2  (Concluded) 


Salinity,  ppt 


Trinity 

Hydrology  Months  Point  2 


Channel 

East  Bay 

South 

West 

West  Bay 

Point  3 

Point  7 

Point  8 

Point  12 

2049  High 


2024  Low 


2024  High 


2049  Med 


2049  High 


Apr-May  0 

Jun-Jul  0 

Aug-Sep  0 

Apr-May  0 

Jun-Jul  0 

Aug-Sep  0 

Apr-May  0 

Jun-Jul  1 

Aug-Sep  1 

Apr-May  1 

Jun-Jul  2 

Aug-Sep  2 

Apr-May  0 

Jun-Jul  0 

Aug-Sep  1 

Apr-May  0 

Jun-Jul  0 

Aug-Sep  1 

Apr-May  0 

Jun-Jul  1 

Aug-Sep  1 

Apr-May  0 

Jun-Jul  0 

Aug-Sep  1 

Apr-May  0 

Jun-Ju!  0 

Aug-Sep  1 


Phase  1  (Continued) 


Phase  II 
1 
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the  filtering  process  the  first  and  last  28-day  periods  show  an  artificial  bulge, 
and  therefore  these  were  not  included  in  Table  2. 

The  typical  Galveston  Bay  trend  is  for  low  salinity  values  through  early 
summer.  The  large  drop  in  freshwater  inflow  beginning  in  July  along  with  the 
higher  Gulf  salinity  causes  a  rise  in  salinity,  usually  reaching  a  maximum 
around  October.  The  model  Gulf  salinity  drops  quickly  in  September  and  the 
bay  salinities  follow  later.  The  response  to  the  freshwater  changes  are,  of 
course,  more  pronounced  for  the  high-fiow  run  than  at  lower  flow  conditions. 

The  overall  increase  in  salinity  of  the  Phase  I  deepening  in  comparison  to 
the  existing  conditions  was  regionally  most  noticeable  in  the  upper  bay  and  in 
the  fall  for  all  regions.  During  high  freshwater  flow  periods,  this  region  and 
Trinity  Bay  were  practically  fresh  for  both  existing  and  Phase  I.  Generally, 
the  maximum  increase  in  salinity  was  greater  during  the  high-flow  scenario 
than  in  the  low.  The  salinity  gradient  from  the  channel  to  the  bay  surrounding 
it  was  larger  as  a  result  of  the  deepening.  The  increase  in  the  lower  bay  was 
much  smaller  than  in  other  areas,  typically  about  1  ppt. 

Considering  first  the  medium-flow  conditions,  the  Trinity  Bay  area 
(points  1,  2,  13  and  14)  shows  a  small  increa.se  in  maximum  salinity  in  the  late 
summer  and  fall  on  the  order  of  1  or  2  ppt.  During  the  high-flow  period  in 
late  spring  and  early  summer  the  difference  was  less.  The  upper  bay  (points  3, 
4  and  15)  shows  the  most  dramatic  increase;  in  fact,  point  15  has  only  about 
2  months  when  the  salinity  was  less  than  10  ppt  versus  a  6-month  period  for 
the  existing  conditions.  The  maximum  increase  in  salinity  for  this  location 
was  on  the  order  of  5  ppt  occurring  in  the  fall.  The  largest  change  for  the 
midbay  (locations  6,  9,  and  10)  was  an  increase  of  about  2  ppt  in  the  fall. 

East  and  West  Bays  (locations  7  and  12,  respectively)  are  only  slightly  affected 
by  the  deepening  with  West  Bay  showing  about  a  1-ppt  increase  in  late  sum¬ 
mer  and  fall.  The  lower  bay  (location  8)  in  the  early  summer  actually  shows  a 
slight  decrease  in  salinity  due  to  Phase  I.  The  channel  is  somewhat  more 
stratified  with  the  enlarged  channel  perhaps  leading  to  less  saline  flows  leaving 
the  channel  in  the  lower  bay. 

The  low  flow  shows  a  more  uniform  plan  impact  over  the  year,  and  con¬ 
versely  the  high  flow  shows  a  less  uniform  change.  The  upper  bay  does 
slightly  exceed  20  ppt  for  Phase  I,  whereas  the  base  existing  conditions  are 
about  17  or  18  ppt  throughout  the  fall.  The  increase  due  to  the  plan  is  gener¬ 
ally  greater  in  the  high-flow  year  than  in  low;  thus,  the  range  of  variation  from 
low  to  high  is  less  with  the  Phase  I  condition.  These  are  shown  in  Plates  49- 
160,  respectively.  The  salinity  time-series  for  NED  and  Phase  II  are  shown  in 
Plates  161-176  and  177-272. 

These  future  year  hydrology  scenarios  for  2024  and  2049  were  tested  for  a 
low-,  medium-,  and  high-flow  condition  and  for  9  months  (January- 
September).  These  tests  were  conducted  for  existing  and  for  Phase  I  and  II 
bathymetry/geometry.  The  future  year  hydrology  indicated  as  1999  was  tested 
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for  the  same  9-monlh  period  for  low-flow  conditions  and  existing  and  Phase  I 
bathymetry/geometry. 

The  change  in  salinity  between  Phase  I  and  existing  bathymetry/geometry 
appears  to  be  less  significant  in  all  future  hydrologic  scenarios  (1999,  2024, 
and  2049).  In  all  these  hydrologies,  the  flow  distribution  increases  flow  from 
Buffalo  Bayou  and  the  San  Jacinto  Rivers  and  reduces  flow  from  Trinity 
River.  So  the  Trinity  Bay  area  has  increased  salinity  under  existing 
bathymetry.  The  additional  flow  to  the  west  side  of  the  bay  produces  a 
decrease  in  salinity  there.  The  deepened  channel  results  in  less  increase  in 
salinity  in  the  future  hydrologies. 

The  salinity  increases  are  greatest  for  the  2024  Phase  II  scenario.  These  are 
larger  than  the  1990  Phase  I  scenario  reported  previously.  In  2024  the 
distribution  of  change  is  similar  to  that  reported  previously,  though  the 
difference  is  noticeable  throughout  the  spring  and  summer.  The  overall 
differences  for  the  year  2049  were  much  less;  in  fact,  they  appear  to  be  less 
than  for  the  Phase  I  conditions  under  1990  hydrology.  In  the  lower  bay  it  is 
interesting  that  for  a  period  of  time  in  the  spring  the  existing  condition  has  a 
higher  salinity  than  Phase  11.  This  has  been  noted  in  earlier  studies  (Berger 
and  Boland  1979)  in  bay  regions  away  from  the  channel.  It  is  likely  that  the 
deepened  and  more  flood  dominant  channel  allows  less  fresh  water  to  exit 
Galveston  Bay  through  the  channel  and  thus  forces  the  fresh  water  through  the 
bay.  The  channel  is  also  slightly  more  stratified  under  deepened  conditions. 
During  the  period  of  salinity  rebound,  the  salinity  in  the  lower  bay  appears  to 
react  more  slowly,  perhaps  due  to  the  presence  of  more  stratification.  It  may 
be  that  more  salt  is  carried  in  the  channel  past  the  lower  bay  rather  than 
immediately  raising  lower  bay  salinity. 


Isohaline  plots 

The  data  were  further  presented  as  isohaline  (lines  of  equal  salinity)  plots. 
In  Plates  273-308,  the  lines  depict  equal  salinities  for  the  existing  versus 
Phase  I  bathymetries  in  5-ppt  intervals.  The  low-flow  case  is  shown  in 
Plates  273-284,  medium  flow  in  285-296,  and  high  flow  in  297-308.  Gener¬ 
ally  speaking,  the  deepened  channel  in  the  Phase  I  bathymetry  shifts  the  isoha¬ 
lines  further  up  the  channel  and  deeper  into  Trinity  Bay,  with  the 
August-October  shifts  being  the  most  pronounced.  The  medium-  and  high- 
flow  cases  indicate  similar  patterns  for  both  bathymetries  as  depicted  for  the 
low  flow  for  the  months  of  January  through  September.  Larger  differences  are 
observed  for  October  through  December. 

The  existing  versus  Phase  I  case  for  1999  hydrology  is  shown  in 
Plates  309-317;  this  is  a  low-flow  scenario.  The  existing  versus  Phase  I  cases 
are  shown  for  low  (Plates  318-326),  medium  (Plates  327-335),  and  high 
(Plates  336-344)  fresh  water,  2024  hydrology.  These  show  a  comparison 
between  plan  and  existing  that  is  much  closer  than  for  the  1990  hydrology. 

For  the  low-  and  medium-inflow  conditions  the  most  notable  differences  are  in 
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the  upper  west  side  of  the  bay  beginning  in  May  and  continuing  through 
September.  For  the  high-flow  condition  these  differences  begin  in  August. 

The  existing  versus  Phase  I  cases  are  shown  for  low  (Plates  345-353), 
medium  (Plates  354-362),  and  high  (Plates  363-371)  fresh  water, 

2049  hydrology.  These  results  are  quite  similar  to  the  2024  hydrology. 

The  existing  versus  Phase  II  cases  are  shown  for  low  (Plates  372-380), 
medium  (Plates  381-389),  and  high  (Plates  390-398)  freshwater  inflow, 

2024  hydrology.  These  depict  slightly  increased  salinity  intrusion  for  the 
months  January  through  May,  and  larger  increases  for  the  drier  months  of  June 
through  September.  The  pattern  is  similar  to  the  existing  hydrology.  Phase  I 
comparison  to  the  base  condition.  This  low-flow  hydrology  also  represents  the 
lowest  of  the  low-flow  hydrologies  tested,  and  thus  should  represent  the  “worst 
case”  as  far  as  salinity  intrusion  results,  as  comparison  to  Plates  318-344 
indicates. 

The  existing  versus  Phase  II  cases  are  shown  for  low  (Plates  399-407), 
medium  (Plates  408-416),  and  high  (Plates  417-425)  freshwater  inflow, 

2049  hydrology.  These  depict  improved  conditions  when  compared  to  the 
2024  cases,  and  are  very  similar  to  the  results  for  Phase  I,  1990  hydrology. 
This  can  be  seen  by  comparing  Plates  399-425  with  Plates  345-371. 


Mitigation  plan 

The  channel  configuration  for  this  test  is  identical  to  the  Phase  I  project, 
except  for  the  addition  of  the  Bayport  Channel  and  a  small  channel  through 
Atkinson  Island,  and  the  final  Beneficial  Uses  Group  dredged  material  place¬ 
ment  plan  was  used  instead  of  the  preliminary  plan,  which  consisted  of  the 
17  islands.  The  conditions  under  which  this  bathymetry  was  tested  were  for 
the  2024  year  hydrology  distribution  for  a  typical  medium  freshwater  inflow 
condition.  This  run  was  for  1  year.  As  expected,  these  results  are  very  much 
like  the  Phase  I  results  for  the  same  hydrology  and  inflow;  therefore,  the  miti¬ 
gation  plan  will  have  no  greater  impact  on  salinities  than  Phase  I. 
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4  Summary  and  Conclusions 


In  this  report  the  results  of  extensive  model  testing  of  a  proposed  enlarge¬ 
ment  of  the  Houston-Galvcslon  Navigation  Channels  through  Galveston  Bay, 
Texas,  are  presented.  The  present  channel  project  depth  is  40  ft  and  the  deep¬ 
ened  project  depths  are  45  ft  (Phase  1)  and  50  ft  (Phase  II).  Several  How 
distributions  within  the  basin  were  tested.  These  arc  referred  to  as  the  1990, 
1999,  2024,  and  2049  year  hydrologies.  Within  each  of  these,  a  typical  low-, 
medium-,  and  high-llow  year  was  developed.  The  future  hydrologic  years 
represent  predicted  redistribution  or  increase  of  freshwater  tlow  within  the 
basin.  The  1999  year  is  a  result  of  How  redistribution  that  might  be  expected 
from  implementation  of  the  Wallisville  Dam.  The  years  2024  and  2049  repre¬ 
sent  the  increased  water  demand  of  the  city  of  Houston,  the  results  of  which 
are  for  more  fresh  water  to  enter  the  bay  via  Buffalo  Bayou  and  the  San 
Jacinto  River. 

Results  of  testing  for  the  45-ft  channel  showed  that  largest  increases  in 
salinity  were  in  low-salinity  areas,  including  the  upper  west  side  of  the  bay 
across  the  channel  from  Atkinson  Island,  and  the  upper  bay  channel.  Trinity 
Bay  showed  a  small  salinity  increase.  South  of  a  line  from  Smith  Point  to 
Eagle  Point  the  increase  was  generally  less  than  1  ppt.  Some  locations  in  the 
south  bay  near  the  navigation  channel  occasionally  showed  a  decrease  in 
salinity  for  the  deeper  channel  configuration.  These  decreases  occurred  during 
the  period  of  rebound  in  salinity  after  the  high  inllow  period  of  late  spring. 

The  50-ft  channel  resulLs  were  similar,  qualitatively,  to  those  of  the  45-ft 
channel.  However,  the  increases  in  salinity  relative  to  existing  channel 
conditions  were  larger. 

The  deepened  channels  showed  increased  salinity  stratification.  The  strati¬ 
fication  increased  with  channel  project  depths  and  with  freshwater  inflow  in 
the  Buffalo  Bayou/San  Jacinto  River  Basin.  The  slight  decrease  in  salinity 
under  deepened  channel  conditions  in  the  southern  part  of  the  bay  is  probably 
a  result  of  this  increased  stratification.  The  higher  salinity  may  be  propagating 
up  the  channel  rather  than  spilling  out  into  the  adjacent  lower  bay. 

The  1999,  2024,  and  2049  hydrologic  years  result  in  more  freshwater 
inflow  to  Galveston  Bay  through  the  Buffalo  Bayou  and  San  Jacinto  River, 

The  deepened  channel  typically  resulted  in  less  significant  salinity  increases  in 
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these  future  scenarios  than  for  the  1990  year.  These  future  scenarios  redistri¬ 
bute  some  of  the  freshwater  inflow  from  the  Trinity  River  and  reintroduce  it 
through  Buffalo  Bayou  and  San  Jacinto  River.  The  model  shows  a  cor¬ 
responding  increase  in  salinity  in  Trinity  Bay  and  decrease  in  the  wester  upper 
bay  salinity. 

The  mitigation  plan  was  the  same  channel  geometry  as  that  of  Phase  I  but 
with  fewer  and  a  different  arrangement  of  disposal  islands.  The  tests  indicated 
that  the  results  were  quite  similar  to  those  of  Phase  I. 
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FEBRUARY  1990,  PHASE  I 
LOW  FLOW,  BOTTOM 


San  Jadnu)  River 


SALINITY  CONTOURS,  PPT 
MARCH  1990,  PHASE  I 
LOW  FLOW.  BOTTOM 


Plate  275 


SALINITY  CONTOURS,  PPT 
APRIL  1990,  PHASE  I 
LOW  FLOW,  BOTTOM 

Plate  276 


San  Jacimo  River 


Triniiv  River 


SALINITY  CONTOURS,  PPT 
MAY  1990,  PHASE  I 
LOW  FLOW,  BOTTOM 


Plate  277 


Plate  278 


SALINITY  CONTOURS,  PPT 
JUNE  1990,  PHASE  I 
LOW  FLOW,  BOTTOM 


San  Jacinu)  River 


Bulfalo  Bayou 


I  riniiv  River 


Clear  Lake^ 


Smith  Point 


Eaaie  Point 


Texas  City 


Rollover 

Pass 


Base 
Phase  I 


Galveston  f 


San  Luis  Pass 


SALINITY  CONTOURS,  PPT 
JULY  1990,  PHASE  I 
LOW  FLOW,  BOTTOM 


Plate  279 


SALINITY  CONTOURS,  PPT 

AUGUST  1990,  PHASE  I 
LOW  FLOW,  BOTTOM 


SALINITY  CONTOURS,  PPT 
SEPTEMBER  1990,  PHASE  I 
LOW  FLOW,  BOTTOM 


San  Jaciiuo  River 


San  Luis  Puss 


SALINITY  CONTOURS,  PPT 
NOVEMBER  1990,  PHASE  1 
LOW  FLOW,  BOTTOM 


San  Jadnu)  River 


SALINITY  CONTOURS.  PPT 
FEBRUARY  1990,  PHASE  I 
MEDIUM  FLOW,  BOTTOM 


San  JacinU)  River 


SAUNITY  CONTOURS,  PPT 
APRIL  1990,  PHASE! 
MEDIUM  FLOW.  BOTTOM 


Plate  288 


SALINITY  CONTOURS,  PPT 
MAY  1990,  PHASE  1 
MEDIUM  FLOW,  BOTTOM 


SAUNITY  CONTOURS,  PPT 
JULY  1990,  PHASE  I 
MEDIUM  FLOW,  BOTTOM 


Plate  296 


SALINITY  CONTOURS.  PPT 
DECEMBER  1990,  PHASE  I 
MEDIUM  FLOW,  BOTTOM 


SALINITY  CONTOURS,  PPT 
FEBRUARY  1990,  PHASE  I 
HIGH  FLOW,  BOTTOM 


SALINITY  CONTOURS,  PPT 
MARCH  1990,  PHASE  I 
HIGH  FLOW,  BOTTOM 


SALINITY  CONTOURS,  PPT 
APRIL  1990,  PHASE  I 
HIGH  FLOW.  BOTTOM 


SALINITY  CONTOURS,  PPT 
JUNE  1990,  PHASE  I 
HIGH  FLOW,  BOTTOM 


SALINITY  CONTOURS,  PPT 
JULY  1990,  PHASE  I 
HIGH  FLOW,  BOTTOM 


San  Jacimo  River 


Plate  304 


SALINITY  CONTOURS,  PPT 
AUGUST  1990,  PHASE  I 
HIGH  FLOW,  BOTTOM 


SALINITY  CONTOURS,  PPT 
OCTOBER  1990.  PHASE  I 
HIGH  FLOW,  BOTTOM 


Plate  307 


SALINITY  CONTOURS,  PPT 
JANUARY  1999,  PHASE  I 
LOW  FLOW,  BOTTOM 


Plate  309 


San  Jacinto  River 


Plate  310 


SALINITY  CONTOURS,  PPT 
FEBRUARY  1999,  PHASE! 
LOW  FLOW.  BOTTOM 


SALINITY  CONTOURS,  PPT 
MARCH  1999,  PHASE  I 
LOW  FLOW,  BOTTOM 


Plate  311 


San  Jacinto  River 


Plate  312 


SALINITY  CONTOURS,  PPT 
APRIL  1999,  PHASE  I 
LOW  FLOW,  BOTTOM 


San  Jacinio  River 


Buffalo  Bavou 


Trinity  River 


Clear  Lake 


Smith  Point 


Eagle  Point 


Texas  City 


Rollover 

Pass 


Base 
Phase  I 


Galveston) 


San  Luis  Pass 


SALINITY  CONTOURS,  PPT 
MAY  1999,  PHASE  I 
LOW  FLOW,  BOTTOM 


Plate  313 


San  Jacinto  River 


Buffalo  Bavou 


Trinity  River 


N 


Clear  Lake^ 


Smith  Point 


Eagle  Point 


20^  20" 


Rollover 

Pass 


Texas  City 


Base 
Phase  I 


Galveston  f 


San  Luis  Pass 


SALINITY  CONTOURS,  PPT 
JULY  1999.  PHASE  I 
LOW  FLOW,  BOTTOM 


Plate  315 


SALINITY  CONTOURS,  PPT 
AUGUST  1999,  PHASE  I 
LOW  FLOW,  BOTTOM 


San  Jacinto  River 


SALINITY  CONTOURS.  PPT 
SEPTEMBER  1999,  PHASE  I 
LOW  FLOW,  BOTTOM 


Plate  317 


Plate  318 


SALINITY  CONTOURS,  PPT 
JANUARY  2024,  PHASE  I 
LOW  FLOW.  BOTTOM 


San  Jacinto  River 


Buffalo  Bavou 


Trinity  River 


Clear  Lake 


Smith  Point 


Eagle  Point 


Texas  City 


Rollover 

Pass 


Base 
Phase  I 


Galveston! 


San  Luis  Pass 


SALINITY  CONTOURS,  PPT 
FEBRUARY  2024,  PHASE  I 
LOW  FLOW.  BOTTOM 


Plate  319 


SALINITY  CONTOURS,  PPT 
MARCH  2024,  PHASE  I 
LOW  FLOW.  BOTTOM 


San  Jacinto  River 


Plate  322 


SALINITY  CONTOURS,  PPT 
JUNE  2024,  PHASE  I 
LOW  FLOW,  BOTTOM 


Plate  323 


SALINIPi'  CONTOURS,  PPT 
AUGUST  2024,  PHASE  i 
LOW  FLOW,  BOTTOM 


Plate  325 


Plate  326 


SALINITY  CONTOURS,  PPT 
SEPTEMBER  2024,  PHASE  I 
LOW  FLOW,  BOTTOM 


SALINITY  CONTOURS,  PPT 
JANUARY  2024,  PHASE  I 
MEDIUM  FLOW.  BOTTOM 


Plate  327 


SALINITY  CONTOURS.  PPT 
FEBRUARY  2024.  PHASE  I 
MEDIUM  FLOW,  BOTTOM 


Trinity  River 


SALINITY  CONTOURS,  PPT 
MARCH  2024.  PHASE  I 
MEDIUM  FLOW,  BOTTOM 


Plate  329 


SALINITY  CONTOURS,  PPT 
MAY  2024,  PHASE  1 
MEDIUM  FLOW,  BOTTOM 


Plate  331 


San  Jacinto  River 


SALINITY  CONTOURS,  PPT 
JUNE  2024,  PHASE  I 
MEDIUM  FLOW,  BOTTOM 


\ 


SALINITY  CONTOURS,  PPT 
JULY  2024,  PHASE  I 
MEDIUM  FLOW,  BOTTOM 


Plate  333 


Plate  334 


SALINITY  CONTOURS,  PPT 

AUGUST  2024.  PHASE  I 
MEDIUM  FLOW,  BOTTOM 


SALINITY  CONTOURS,  PPT 
SEPTEMBER  2024,  PHASE  1 
MEDIUM  FLOW,  BOTTOM 


Plate  335 


25 


SALINITY  CONTOURS,  PPT 
FEBRUARY  2024,  PHASE  I 
HIGH  FLOW.  BOTTOM 


San  Jacinto  River 


Buffalo  Bavou 


Trinity  River 


Clear  LakCr- 


Eagle  Point^ 


Smith  Point 


Texas  City 


Base 
Phase  I 


Galveston 


Rollover 

Pass 


SALINITY  CONTOURS.  PPT 
MARCH  2024.  PHASE  I 
HIGH  FLOW.  BOTTOM 


SALINITY  CONTOURS.  PPT 
JUNE  2024,  PHASE  I 
HIGH  FLOW.  BOTTOM 


Plate  341 


SALINITY  CONTOURS,  PPT 
AUGUST  2024.  PHASE  I 
HIGH  FLOW.  BOTTOM 


SALINITY  CONTOURS.  PPT 
SEPTEMBER  2024,  PHASE  I 
HIGH  FLOW.  BOTTOM 


SALINITY  CONTOURS.  PPT 
JANUARY  2049,  PHASE  I 
LOW  FLOW,  BOTTOM 


Plate  345 


SALINITY  CONTOURS,  PPT 
MARCH  2049.  PHASE  1 
LOW  FLOW,  BOTTOM 


Plate  347 


SALINITY  CONTOURS,  PPT 

APRIL  2049,  PHASE  I 
LOW  FLOW.  BOTTOM 


San  Jacimo  River 


SALINITY  CONTOURS,  PPT 
MAY  2049,  PHASE  I 
LOW  FLOW,  BOTTOM 


Plate  349 


SALINITY  CONTOURS,  PPT 
JULY  2049,  PHASE  I 
LOW  FLOW,  BOTTOM 


Plate  351 


SALINITY  CONTOURS,  PPT 
AUGUST  2049,  PHASE  I 
LOW  FLOW,  BOTTOM 


San  Jacinto  River 


Buffalo  Bavou 


Trinity  River 


Clear  Lake^ 


Smith  Point 


Eagle  Point 


Texas  City 


Rollover 

Pass 


Base 
Phase  I 


Galveston/ 


San  Luis  Pass 


SALINITY  CONTOURS,  PPT 
SEPTEMBER  2049,  PHASE  ! 
LOW  FLOW,  BOTTOM 


San  Jacinto  River 


SALINITY  CONTOURS,  PPT 
JANUARY  2049,  PHASE  I 
MEDIUM  aOW.  BOTTOM 


Plate  354 


San  Jacinto  River 


SALINITY  CONTOURS,  PPT 
FEBRUARY  2049.  PHASE  I 
MEDIUM  FLOW.  BOTTOM 


Plate  355 


SAUNITY  CONTOURS,  PPT 
APRIL  2049  PHASE  I 
MEDIUM  FLOW,  BOTTOM 


Plate  357 


San  Jacinto  River 


Plate  358 


SALINITY  CONTOURS,  PPT 
MAY  2049,  PHASE  I 
MEDIUM  FLOW,  BOTTOM 


Plate  359 


SAUNITY  CONTOURS,  PPT 
JULY  2049,  PHASE  I 
MEDIUM  FLOW,  BOTTOM 


SALINITY  CONTOURS.  PPT 
AUGUST  204a,  PHASE  I 
MEDIUM  FLOW,  BOTTOM 


Plate  361 


SALINITY  CONTOURS,  PPT 
SEPTEMBER  2049,  PHASE  I 
MEDIUM  FLOW,  BOTTOM 


Plate  366 


SALINITY  CONTOURS,  PPT 

APRIL  2049,  PHASE  I 
HIGH  FLOW,  BOTTOM 


SALINITY  CONTOURS,  PPT 
MAY  2049,  PHASE  I 
HIGH  FLOW,  BOTTOM 


Plate  367 


San  Jacinto  River 


SALINITY  CONTOURS,  PPT 
JUNE  2049,  PHASE  I 
HIGH  FLOW,  BOTTOM 

Plate  368  - - 


SALINITY  CONTOURS,  PPT 
AUGUST  2049,  PHASE  I 
HIGH  FLOW,  BOTTOM 


San  Jacinto  River 


Trinity  River 


SALINITY  CONTOURS,  PPT 
SEPTEMBER  2049,  PHASE  I 
HIGH  FLOW,  BOTTOM 


SALINITY  CONTOURS.  PPT 
APRIL  2024,  PHASE  II 
LOW  FLOW,  BOTTOM 


SALINITY  CONTOURS,  PPT 
JUNE  2024,  PHASE  II 
LOW  FLOW,  BOTTOM 


Plate  377 


SAUNITY  CONTOURS,  PPT 
AUGUST  2024,  PHASE  II 
LOW  FLOW,  BOTTOM 


San  Jacinto  River 


SALINITY  CONTOURS,  PPT 
JANUARY  2024,  PHASE  II 
MEDIUM  FLOW,  BOTTOM 


Plate  381 


San  Jacinto  River 


SALINITY  CONTOURS,  PPT 
MARCH  2024,  PHASE  II 
MEDIUM  FLOW,  BOTTOM 


Plate  384 


SAUNITY  CONTOURS,  PPT 
APRIL  2024,  PHASE  II 
MEDIUM  FLOW,  BOTTOM 


San  Jacinto  River 


Trinity  River 


SALINITY  CONTOURS,  PPT 
MAY  2024,  PHASE  II 
MEDIUM  FLOW,  BOTTOM 


Plate  385 


SALINITY  CONTOURS,  PPT 
JUNE  2024,  PHASE  II 
MEDIUM  FLOW,  BOTTOM 


San  Jacinto  River 


Trinity  River 


SALINITY  CONTOURS,  PPT 
JULY  2024,  PHASE  II 
MEDIUM  FLOW,  BOTTOM 


Plate  387 


San  Jacinto  River 


Plate  388 


SALINITY  CONTOURS,  PPT 

AUGUST  2024,  PHASE  II 
MEDIUM  FLOW,  BOTTOM 


SALINITY  CONTOURS,  PPT 
SEPTEMBER  2024,  PHASE  II 
MEDIUM  FLOW ,  BOTTOM 


SALINITY  CONTOURS,  PPT 
FEBRUARY  2024,  PHASE  II 
HIGH  FLOW,  BOTTOM 


Plate  391 


Plate  392 


SALINITY  CONTOURS,  PPT 
MARCH  2024,  PHASE  II 
HIGH  FLOW,  BOTTOM 


San  Jacinto  River 


Trinity  River 


SALINITY  CONTOURS,  PPT 
APRIL  2024,  PHASE  II 
HIGH  FLOW,  BOTTOM 


Plate  393 


Plate  394 


SALINITY  CONTOURS,  PPT 

MAY  2024,  PHASE  II 
HIGH  ROW,  BOTTOM 


San  Jacinto  River 


Trinity  River 


Plate  395 


San  Jacinto  River 


Buffalo  Bayou 


Trinity  River 
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Clear  Lake 


Smith  Point 


Eagle  Point 


Texas  City  2C 


Rollover 

Pass 


Base 
Phase  II 


Galveston) 


f20^  3C 


San  Luis  Pass 


SALINITY  CONTOURS,  PPT 
AUGUST  2024,  PHASE  II 
HIGH  FLOW,  BOTTOM 


Plate  397 


San  Jacinto  River 


Plate  398 


SALINITY  CONTOURS,  PPT 
SEPTEMBER  2024,  PHASE  II 
HIGH  FLOW,  BOTTOM 


SALINITY  CONTOURS,  PPT 
JANUARY  2049,  PHASE  II 
LOW  FLOW,  BOTTOM 


SAUNITY  CONTOURS,  PPT 
MARCH  2049,  PHASE  II 
LOW  FLOW,  BOTTOM 


Triniiv  River 

■ 


SALINITY  CONTOURS,  PPT 

JUNE  2049,  PHASE  II 
LOW  FLOW,  BOTTOM 


SALINITY  CONTOURS,  PPT 
JUNE  2049,  PHASE  II 
LOW  FLOW.  BOTTOM 


San  JaciiiU)  River 


SALINITY  CONTOURS,  PPT 
MAY  2049,  PHASE  II 
LOW  FLOW,  BOTTOM 


Plate  403 


Plate  406 


SALINITY  CONTOURS,  PPT 
AUGUST  2049,  PHASE  II 
LOW  FLOW,  BOTTOM 


San  Jacinto  River 


SALINITY  CONTOURS,  PPT 
JANUARY  2049,  PHASE  II 
MEDIUM  FLOW,  BOTTOM 

Plate  408 


San  Jacinto  River 


SALINITY  CONTOURS,  PPT 
FEBRUARY  2049,  PHASE  II 
MEDIUM  FLOW,  BOTTOM 


Plate  409 


San  Jacinto  River 


SALINITY  CONTOURS,  PPT 
MARCH  2049,  PHASE  II 
MEDIUM  FLOW,  BOTTOM 


San  Jacinto  River 


Plate  412 


SALINITY  CONTOURS,  PPT 
MAY  2049,  PHASE  II 
MEDIUM  FLOW,  BOTTOM 


San  Jacinto  River 


Trinity  River 


SALINITY  CONTOURS,  PPT 
JUNE  2049,  PHASE  II 
MEDIUM  FLOW,  BOTTOM 


Plate  413 


Plate  414 


SALINITY  CONTOURS,  PPT 
JULY  2049,  PHASE  II 
MEDIUM  FLOW,  BOTTOM 


SALINITY  CONTOURS,  PPT 
FEBRUARY  2049,  PHASE  II 
HIGH  FLOW,  BOTTOM 


Plate  420 


SALINITY  CONTOURS,  PPT 
APRIL  2049,  PHASE  II 
HIGH  FLOW,  BOTTOM 


San  Jacinto  River 


SALINITY  CONTOURS,  PPT 
MAY  2049,  PHASE  II 
HIGH  FLOW,  BOTTOM 


Plate  421 


Plate  422 


SALINITY  CONTOURS,  PPT 
JUNE  2049,  PHASE  II 
HIGH  FLOW,  BOTTOM 


SALINITY  CONTOURS,  PPT 
AUGUST  2049,  PHASE  II 
HIGH  FLOW,  BOTTOM 


San  Jacinto  River 


SALINITY  CONTOURS,  PPT 
SEPTEMBER  2049,  PHASE  II 
HIGH  FLOW,  BOTTOM 


Plate  425 


Appendix  A 

The  Hydrodynamic  Code 


The  geometric  complexity  of  this  estuary,  with  its  navigation  channel,  mul¬ 
tiple  inlets,  and  many  proposed  disposal  islands,  requires  a  numerical  model 
that  relies  upon  an  unstructured  computational  mesh.  The  code  chosen  is  the 
Galerkin-based  finite  element  model  RMAIO-WES,  which  is  a  U.S.  Army 
Engineer  Waterways  Experiment  Station  (WES)  adaptation  of  the  RMA-10 
code  developed  by  King  (1993).  This  code  computes  time-varying  open- 
channel  flow  and  salinityAemperature  transport  in  1,  2,  and  3  dimensions.  It 
invokes  the  hydrostatic  pressure  and  mild  slope  assumption.  Vertical  turbu¬ 
lence  is  supplied  using  a  Mellor-Yamada  Level  II  (Mcllor  and  Yamada  1982) 
approach  modified  for  stratification  by  the  method  of  Henderson-Sellers 
(1984).  The  salinity/density  relationship  is  based  upon  Pritchard  (1982). 

The  full  three-dimensional  equations  are  reduced  to  a  set  of  two  momentum 
equations,  an  integrated  continuity  equation,  a  convection-diffusion  equation, 
and  an  equation  of  state.  The  simplification  is  a  result  of  the  hydrostatic 
pressure  approximation. 
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Elevation-related  terms  are  defined  in  Figure  Al. 
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Figure  A1 .  Definitions  for  elevation  terms 
where 
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p  =  density 

=  x,y,z  velocity  components 
t  =  time 
P  -  pressure 


r,.  =  pQv 


/  2  2^(l/2) 

+  ^’a) 


+  ip  \V~  cos  (0) 


r  =  -pQ.u  - 


/  2  2n(1/2) 

9gv„{u,  +  v„) 


+  T(.'  W-  sin  (0) 


Q  =  2ojsin((l)) 

CO  =  rate  of  angular  rotation  of  the  earth 
([)  =  local  latitude 
g  =  gravitational  acceleration 
C  =  Chezy  or  Manning  friction  formulation 
i|.'  =  a  coefficient  from  Wu  (1980) 

W  =  wind  speed 

0  =  wind  direction  counterclockwise  from  easterly 
h  =  depth 

=  x,y  velocity  components  at  the  water  surface 
C  =  water  surface  elevation 
=  x,y  velocity 
a  =  bed  elevation 


A2 
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5  =  salinity 

=  diffusion  coefficient  for  salt 
E  =  eddy  viscosity  components 

The  continuity  equation 


du  =  0 

lx  ly 


is  solved  as  a  second  part  of  each  solution  step.  Equation  A6  is  converted  to 
an  appropriate  boundary  value  problem  through  differentiation  with  respect  to 
z.  After  rearrangement  it  takes  the  form 


d  / ^ 
dz  dx  dy 


subject  to  boundary  conditions  specified  for  the  water  surface  and  the  bed. 

w  =  Uy  —  +  Vy  —  +  —  at  the  water  suface 

^  ^  dx  ^  dy  dt 


da  da 

w  -  u  —  +  — 

^  “  dx  “  dy 


at  the  bed 


Note  that  in  these  equations  the  values  of  ii  and  v  will  be  known  at  all  loca¬ 
tions  from  the  previous  part  of  the  solution  step.  Values  of  w  in  this  solution 
are  used  in  the  next  iteration  for  u,v,h,  and  s. 

The  geometric  system  varies  with  time;  i.e.,  the  water  depth  h  varies  during 
the  simulation.  In  order  to  develop  an  Eulerian  form  for  the  solution,  it  is 
desirable  to  transform  this  system  to  one  that  can  be  described  with  a  constant 
geometric  structure.  Early  development  of  the  model  (King  1982)  used  a  o- 
transformation  in  which  the  bed  and  the  water  surface  are  transformed  to  con¬ 
stants.  In  a  later  analysis  of  this  method,  King  (1985)  pointed  out  that  at  loca¬ 
tions  where  a  sharp  break  in  bottom  profile  occurs,  the  transformation  is  not 
unique  and  momentum  in  the  component  directions  may  not  be  correctly  pre¬ 
served.  An  alternative  transformation  that  preserves  the  bottom  profile  as 
defined,  but  transforms  the  water  surface  to  a  constant  elevation  is  now  used 
(z^  transformation). 


This  transformation  is  defined  by: 


x^  =  X 


(AlO) 
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(All) 


(A12) 


where  b  is  the  fixed  vertical  location  to  which  the  water  surface  will  be  trans¬ 
formed.  Equations  A1-A6  and  A7-A9  then  incorporate  the  transformation 
(A10-A12). 

Another  advantage  of  this  transformation  is  that  it  produces  constant 
lines  that  are  close  to  horizontal,  i.e.,  z  =  constant  lines.  This  results  in  less 
fictitious  density-driven  currents  near  bed  profile  breaks  (Stelling  and  van 
Kester  1993).  Since  stratification-related  phenomena  are  usually  nearly  hori¬ 
zontal,  it  is  important  that  the  transformation  leave  constant  surfaces  that  are 
nearly  horizonal.  Considering  the  pressure  gradient  (due  to  the  density 
gradient)  in  this  transformation  produces 


dp 

dP 

dP 

dz^ 

-  = 

-  +  _ 

dx 

dx^ 

clz’ 

dx 

(A13) 


In  a  strongly  stratified  stagnant  system  this  pressure  gradient  should  be  zero. 
However,  note  that  Equation  A13  in  the  transformed  system  is  dependent  upon 
two  terms  (each  of  which  could  be  large)  to  cancel  each  other.  This  could 
cause  artificial  currents  due  to  truncation  and  roundoff  error.  A  transformation 
in  which  clzVclx  =  0,  i.e.,  z''==  z,  will  reduce  this  problem.  Figure  A2  shows 
an  example  for  a  case  similar  to  the  Galveston  project  in  which  a  40-ft-deep 
channel  passes  through  an  8-ft-deep  bay.  Here  b  is  chosen  to  be  an  elevation 
of  0  and  ^  is  2  ft.  Near  the  break  in  the  bed  profile  dz^ldx  is  fairly  small,  or 
z  surfaces  are  nearly  horizontal.  Contrast  this  with  the  o  transformation  in 
Figure  A3.  The  a  =  constant  surfaces  are  far  from  horizontal  along  the  chan¬ 
nel  side  slopes.  The  truncation  and  roundofl  errors  tend  to  drive  fictitious 
currents  that  cau.se  the  denser  salt  water  to  leave  the  channel.  The  z'’ 
transformation  results  in 

-^  =  0(i;  -  J^)  (A14) 


whereas  the  a  transformation  is 

=  o(/0 

Ox 


(A15) 


A4 


which  is  much  larger. 

The  Galerkin  finite  element  approximation  of  Equations  A1-A4  and  A7 
uses  a  quadratic  approximation  for  u,v,w,  and  j  and  linear  for  li  and  P.  The 
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nonlinearity  is  addressed  by  Newton-Raphson  iteration  at  each  time-step. 
Generally  the  iteration  process  is  split  into  calculation  of  Equations  A1-A3, 
then  A7,  followed  by  A4.  This  sequence  is  repeated  until  sufficient 
convergence  is  reached. 
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Table  B1 

Present  Hydrology 


(Continued) 


Table  B1  (Concluded) 


Table  B2 
1999  Hydrology 


11,611.3  11,108.2  13,143.1 


Table  B2  {Concluded) 


Table  B3 
2024  Hydrology 
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2049  Hydrology 


Table  B4  (Concluded) 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  Including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining 
the  data  needed,  and  completing  and  reviewing  the  collectioo  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions 
for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway.  Suite  1204,  Arlington,  VA22202-4302,  arxl  to  the 
Office  of  Marragement  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC20503. _ _ _ 


1 .  AGENCY  USE  ONLY  (Leave  blank) 


REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

September  1995  Report  4  of  a  series 


4.  TITLE  AND  SUBTITLE 

Houston-Galveston  Navigation  Channels,  Texas  Project;  Report  4, 
Three-Dimensional  Numerical  Modeling  of  Hydrodynamics  and  Salinity 

6.  AUTHOR(S) 

R.  C.  Berger,  Robert  T.  McAdory,  Joseph  H.  Schmidt,  William  D.  Martin, 
Larry  H.  Hauck 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Engineer  Waterways  Experiment  Station 
3909  Halls  Ferry  Road,  Vicksburg,  MS  39180-6199 


5.  FUNDING  NUMBERS 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

Technical  Report  HL-92-7 


Jones  and  Neuse  Engineering 
Austin,  TX  77553 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Engineer  District,  Galveston 
Jadwin  Building 
2000  Fort  Point  Road 
I  Galveston,  TX  77550 

11.  SUPPLEMENTARY  NOTES 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


Available  from  National  Technical  Information  Service,  5285  Port  Royal  Road,  Springfield,  VA  22161 . 


12a.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


12b.  DISTRIBUTION  CODE 


1 3,  ABSTRACT  (Maximum  200  words) 

This  report  describes  the  testing  program  conducted  to  evaluate  the  impact  of  enlargement  of  the  Houston-Galveston 
Navigation  Channel  on  the  salinity  and  hydrodynamic  fields  of  tidally  influenced  Galveston  Bay.  The  present  channel 
nominal  dimensions  are  40  ft  deep  at  mean  low  water  (mlw)  and  400  ft  wide.  The  proposed  enlargements  tested  are  for  a 
channel  45  ft  deep  at  mlw  and  530  ft  wide  (Phase  I)  and  50  ft  deep  and  600  ft  wide  (Phase  11).  Current  plans  do  not  include 
the  Phase  II  enlargement.  Salinity  fields  for  these  channel  configurations  and  the  existing  channel  dimensions  are  compared. 
In  a  separate  study  the  results  from  these  simulations  were  used  to  drive  an  ecosystem  model  to  predict  oyster  production. 

Testing  conditions  included  tidal  conditions  and  winds  for  the  year  1984.  The  freshwater  inflows  (developed  outside  this 
study)  were  tested  for  low-,  medium-,  and  high-flow  years.  Additionally,  since  water  demand  in  the  future  is  expected  to 
modify  freshwater  distribution  and  quantities,  future  distributions  for  the  test  year  1999  (Wallisville  Dam  in  place),  2024,  and 


14.  SUBJECT  TERMS 

Estuary 
Galveston  Bay 
Mathematical  model 


Navigation  channels 
Salinity  intrusion 


(Continued) 

15.  NUMBER  OF  PAGES 

482 

16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION  18. 
OF  REPORT 

UNCLASSIFIED 


SECURITY  CLASSIFICATION  19. 
OF  THIS  PAGE 

UNCLASSIFIED 


SECURITY  CLASSIFICATION  20. 
OF  ABSTRACT 


LIMITATION  OF  ABSTRACT 


NSN  7540-01-280-5500 


Standard  Form  298  (Rev.  2-89) 
Prescribed  by  ANSI  Std.  Z39-18 
298-102 


13.  (Concluded) 


The  code  used  (RMAIO-WES)  is  a  Galerkin-based  finite  element  solution  to  simulate  three-dimensional  (3-D) 
unsteady  open-channel  flow.  The  code  represents  3-D  hydrodynamics  using  conservation  of  fluid  mass,  horizontal 
momentum,  and  salinity/temperature  transport  equations  subject  to  the  hydrostatic  assumption. 

Results  of  these  tests  showed  that  the  largest  increases  in  salinity  were  in  low-salinity  areas,  including  the  upper 
west  side  of  the  bay  across  the  channel  from  Atkinson  Island,  and  the  upper  bay  channel.  Trinity  Bay  showed  a 
small  salinity  increase.  South  of  midbay  the  salinity  increases  were  generally  less  than  1  ppt.  Some  locations  in  the 
south  bay  near  the  navigation  channel  occasionally  showed  a  decrease  in  salinity  for  the  deeper  channel 
configurations.  These  decreases  occurred  during  the  period  of  rebound  in  salinity  after  the  high  inflow  period  of  late 
spring. 

The  deepened  channels  showed  increased  salinity  stratification.  The  stratification  increased  with  channel  project 
depths  and  with  freshwater  inflow  in  the  Buffalo  Bayou/San  Jacinto  River  Basin. 

The  future  hydrologic  scenarios  result  in  more  freshwater  inflow  to  Galveston  Bay  through  Buffalo  Bayou  and 
San  Jacinto  River.  The  deepened  channels  typically  resulted  in  less  significant  salinity  increases  in  the  scenarios 
than  for  the  present  hydrologic  year  (1990).  These  future  scenarios  redistribute  some  of  the  freshwater  inflow  from 
the  Trinity  River  and  reintroduce  it  through  Buffalo  Bayou  and  San  Jacinto  River,  The  model  indicates  a 
corresponding  increase  in  salinity  in  Trinity  Bay  and  decrease  in  the  western  upper  bay  salinity. 


